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Abstract General Relativity is able to describe the dynamics of galaxies and larger 
cosmic structures only if most of the matter in the Universe is dark, namely it does 
not emit any electromagnetic radiation. Intriguingly, on the scale of galaxies, there 
is strong observational evidence that the presence of dark matter appears to be nec- 
essary only when the gravitational field inferred from the distribution of the lumi- 
nous matter falls below an acceleration of the order of 10^'" m s^^. In the standard 
model, which combines Newtonian gravity with dark matter, the origin of this ac- 
celeration scale is challenging and remains unsolved. On the contrary, the full set 
of observations can be neatly described, and were partly predicted, by a modifica- 
tion of Newtonian dynamics, dubbed MOND, that does not resort to the existence 
of dark matter On the scale of galaxy clusters and beyond, however, MOND is not 
as successful as on the scale of galaxies, and the existence of some dark matter ap- 
pears unavoidable. A model combining MOND with hot dark matter made of sterile 
neutrinos seems to be able to describe most of the astrophysical phenomenology, 
from the power spectrum of the cosmic microwave background anisotropics to the 
dynamics of dwarf galaxies. Whether there exists a yet unknown covariant theory 
that contains General Relativity and Newtonian gravity in the weak field limit, and 
MOND as the ultra-weak field limit is still an open question. 
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1 Introduction 

The dynamical properties of galaxies and larger cosmic structures should be de- 
scribed by their mass content if gravity is the dominant force field. In 1932, Oort 
noticed a shortage of mass required to describe the velocity of stars in the solar 
neighbourhood 111 1 II . In 1933, Zwicky found the very same problem by applying 
the virial theorem to the Coma cluster, based on the radial velocities of a few galax- 
ies derived from their optical spectra 111531 . These early discoveries were revived in 
the seventies when the existence of flat rotation curves in disk galaxies 11201 and 
the requirement of the dynamical stability of galactic disks 11 121 1441 showed that 
a large fraction of mass, in addition to the observed luminous mass, was necessary 
to describe the dynamics of galaxies. Ever since, the evidence of missing mass in 
cosmic structures has become overwhelming ||39| . 

General Relativity (GR) has encountered such a plethora of successes in the solar 
system that the scientific community finds it hard to suppose that GR, and its New- 
tonian weak field limit, might fail on cosmic scales. It is more natural to suppose 
that this missing mass is actually some form of dark matter (DM) that does not emit 
any electromagnetic radiation but only acts gravitationally. 

This idea meets the demand of particle physics: going beyond the standard model 
of particle physics requires the existence of still unknown elementary particles, like, 
for example, supersymmetric particles, axions, Kaluza-Klein excitations, or sterile 
neutrinos. These particles may naturally play the role of the astrophysical DM 1181 : 
they formed in the early Universe and their relic abundance can fill in the fraction of 
mass that is required to describe both the internal dynamics of structures and their 
formation by gravitational instability. 

In addition, most of these particles have the advantage of being non-baryonic, 
namely, in the astrophysical jargon, they are neither electrons nor particles made of 
quarks. Non-baryonic matter includes neutrinos and the hypothetical Weakly Inter- 
acting Massive Particles (WIMPs). So far, only neutrinos have been detected; all 
other non-baryonic matter is still hypothetical. 

Being non-baryonic is advantageous if gravitational instability drives the forma- 
tion of the cosmic structure. In fact, in the inflationary scenario of the standard hot 
Big Bang cosmology, small perturbations due to quantum fluctuations are inflated to 
cosmic scales by the ~ 100 e-folding expansion of the Universe and provide the ini- 
tial conditions of the matter density field. However, this scenario is contradicted by 
the Cosmic Microwave Background (CMB) anisotropics if the matter density field is 
mostly baryonic. In fact, the observation of temperature anisotropics 5T /T ^ 5 /3, 
with 5 the baryonic matter density fluctuations, in the CMB, which formed by red- 
shift z ^ 10^, yields 8T /T ^ 10^^ on ^ 7° angular scales, corresponding to su- 
perc lusters and larger structures [11351 . Gravitational in stability yields a growth rate 
°= (1 +z)^^ or slower Thus, superclusters would have matter overdensities 8 of 
the order of ^ 10^^ today, rather than the observed 5^1 — 10. Non-baryonic DM 
decouples from the radiation field much earlier than baryons and its density pertur- 
bations can start growing at the time of equivalence, when the radiation and matter 
energy densities are equal. At the time of the baryon-radiation decoupling, DM per- 
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turbations have already grown to 5 ~ 10"^ — 10"^ on supercluster scales and they 
will keep growing to 5 ~ 1 — 10 by today. 

Further evidence of the non-baryonic nature of DM is the abundance of light 
elements which are synthetized in the early Universe. Measures of the primordial 
abundance of deuterium, for example, which is particularly sensitive to the photon- 
to-baryon ratio, implies a baryon density Qbli^ = 0.0214 ± 0.002cQ 169), which 
agrees with the bai-yon density i2fo/r = 0.02229 ± 0.00073 implied by the CMB 
anisotropies II136I . yet it is smaller than the total matter density Q,„h^ = 0. 1 14 fTTll . 

By combining all these pieces of evidence, the standard interpretation of the ob- 
servations of the cosmic structure, from galactic scales to the CMB, suggests that 
only 17% of the matter in the Universe is baryonic, whereas the remaining 83% is 
required to be non-baryonic cold DM |7r|. The entire amount of matter however 
yields only 27% of the matter-energy density required to make the Universe geo- 
metrically flat, as suggested by the CMB power spectrum IIJTI . The missing 73% 
can be easily described by a cosmological constant ||T9l , but a large fraction of the 
cosmological community is searching for more sophisticated models, some of them 
including an additional dark energy (DE) fluid ll39l . 

The current standard ACDM model thus contains non-baryonic cold DM and 
a cosmological constant A . However, we are in the uneasy situation where only 
4.5% of the matter-energy density of the Universe is made of matter we can detect 
in laboratories on Earth. The rest is hypothetical (DM) and has unusual properties, 
like negative pressure (A or DE). To avoid this situation, the alternative solution is to 
assume that GR breaks down at some scale. A number of extended or modified theo- 
ries of gravity have been suggested. Among them, Milgrom ll94ll proposed Modified 
Newtonian Dynamics (MOND), an empirical law that modifies Newtonian dynam- 
ics and actually goes beyond a simple alternative theory of gravity. MOND can ele- 
gantly deal with the galaxy scale phenomenology without DM but is less successful 
at describing the dynamical properties of galaxy clusters. Famaey and McGaugh 
||48ll have recently extensively reviewed this model. Here, we summarize the suc- 
cesses of MOND and mention how the combination of MOND with the existence 
of sterile neutrinos of mass in the range 11 eV — 1 ke V can provide a reasonable 
description of the observed phenomenology, from galaxy scales to the CMB. 



2 The Acceleration Scale 

In 1983, MOND was originally introduced to explain the observed rotation curves 
of spiral galaxies 1941 : the observed velocities of the stars are larger than the ve- 
locities that Newtonian dynamics would predict based on the distribution of the lu- 
minous matter 11201 . As we mentioned above, this observational piece of evidence 
suggested the existence of some hidden mass, or DM, responsible for a stronger 
gravitational field and hence larger velocities. An alternative solution is that, on 



' We use the standard parametrization for the Hubble constant today Hq = 100 h km s ' Mpc 
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these cosmic scales, Newtonian gravity does not hold. Milgrom 1941 went beyond 
the suggestion of a modified theory of gravity that becomes relevant beyond some 
large length scale. He rather suggested that Newtonian dynamics breaks down be- 
low an acceleration scale oq. In the formulation of classical physics, the relation 
between the acceleration a acting on a point mass and the Newtonian acceleration 
aN is 



where jj. (x) is an unknown interpolating function that satisfies the conditions (x — > 
+00) = 1 and — > 0) = jc0 To describe the observed rotation curves of spiral 
galaxies, we need to have the acceleration scale ao ~ 10^'" m s^^ 1951 . Moreover, 
Milgrom showed that this modification of Newtonian dynamics also drastically re- 
duces the amount of DM required in groups and clusters of galaxies ll96l . 

More importantly, Milgrom's suggestion provided clean predictions ll95l that 
were confirmed in later years. Most notably: (1) the zero point of the Tully-Fisher 
relation; (2) the one-to-one correspondence between features in the baryonic distri- 
bution and the rotation curve; (3) a larger mass discrepancy, when interpreted with 
Newtonian gravity, in low surface brightness dwarf spheroidal galaxies; and (4) the 
validity of the Tully-Fisher relation for low surface brightness disk galaxies. 

It also became clear very soon that there are different numerical coincidences that 
might suggest that oq is indeed a fundamental quantity (see fTTl for a recent review). 
For example, it is rather intriguing that oq is related to the Hubble constant Hq and 
the cosmological constant A with the relations aq = cHo/2n and Oq — c^A/2n, 
where c is the speed of light. 

A compilation of kinematic and photometric data of disk galaxies very clearly 
shows the validity of this ansatz on the acceleration scale l87l : disk galaxies do not 
require DM beyond a given length scale, but rather below a given acceleration scale. 
Figure [T] shows the ratio between the observed velocity V and the velocity Vj, ex- 
pected in Newtonian gravity from the distribution of the visible matter The square 
of this ratio is proportional to the ratio between the total mass and the visible mass 
in a spherical system. Evidently, (V/V/,)^ > 1 reflects a mass discrepancy and the 
requirement of DM; the larger the ratio the larger is the contribution of DM. Figure 
[T]clearly shows that the need for DM increases with decreasing Newtonian acceler- 
ation — /r which is proportional to the baryonic surface density. The MOND 
prediction is thus that low-surface brightness (LSB) galaxies are more DM domi- 



^ A year later, Bekenstein and Milgrom (151 showed that, in spherical symmetry, equation ([T) is 
equivalent to a modified theory of gravity where the standard Poisson equation is replaced by 



where 4> is the gravitational potential and p the mass distribution. A more recent modified theory 
of gravity that reproduces equation ([T) is the quasi-linear formulation of MOND (QUMOND) 
11031 . This theory has the advantage of involving only linear differential equations and one non- 
linear algebraic step. See 1481 for a comprehensive review of the various formulations of MOND 
as modified dynamics or modified gravity. 




(1) 
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Fig. 1 (V/Vfc)^ vs. the Newtonian gravity = V(j"/r, derived from the baryonic surface density, 
in almost one hundred spiral galaxies. The need for DM appears when (V/V/,)- > 1. Observations 
clearly show that this only happens below a specific acceleration of the order of 10^'" m s^^. 
Courtesy of B. Famaey and S. McGaugh. Reproduced from 1481 . with permission. 

nated than high-surface brightness (HSB) galaxies in Newtonian gravity. However, 
the Tully-Fisher relation, that we describe below, remains valid for both classes of 
galaxies. In the DM paradigm, it is unclear how the tight correlation shown in Figure 
[T]can emerge from the combination of intrinsically chaotic processes, including the 
history of the DM halo formation, the star formation history and feedback processes. 

If this acceleration scale does indeed exist, we have two relevant consequences: 
(1) MOND has to break the Strong Equivalence Principle, which states that all laws 
of physics are independent of velocity and location in spacetime, and (2) Birkhoff 's 
theorem is not valid. 

MOND is not required to break the Weak Equivalence Principle which asserts 
the equality between inertial and gravitational masses. In GR, the Weak Equiv- 
alence Principle translates into the assumption that all matter fields are coupled 
to the geodesic metric g^y, whereas the Strong Equivalence Principle is guar- 
anteed by assuming that the same metric obeys the Einstein-Hilbert action 
Seh = c'*/(167rG) / g^^ Rij^v^—gY^^d'^x, where R^y is the Ricci tensor and g the 
determinant of the metric. Therefore, one way to break the Strong Equivalence Prin- 
ciple is to keep the Einstein and geodesic metrics distinct. 

The non-validity of Birkhoff's theorem clearly complicates the interpretation of 
the dynamical properties of self-gravitating systems. In Newtonian gravity, any ob- 
ject is subject to the gravitational attraction of all the other objects in the Universe, 
but in Newtonian linear dynamics a star in a galaxy or a galaxy in a galaxy cluster 
are virtually isolated from the gravitational field felt by the parent system, unless this 
latter external field varies across the parent system and originates the well known ef- 
fect of tides. In MOND, the external field affects the internal motions of the objects 
even when this field is constant, because it is the total acceleration felt by the object 
that determines the dynamical, MONDian or Newtonian, regime. We will mention 
below that this issue can be relevant in the investigation of the internal dynamics 
of elliptical galaxies, that are mostly located in regions of high galaxy density, and 
dwarf spheroidal galaxies and globular clusters, that are satellites of larger galaxies. 
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In Sect, m we describe an example of how these issues enter the construction of 
a covariant theory of MOND. 



3 Self-gravitating Systems 
3.1 Disk Galaxies 

On the scales of galaxies there are a number of observations that were either pre- 
dicted or easily explained by MOND. In the DM paradigm, these very same obser- 
vations require extreme fine tuning between baryonic and non-baryonic physics or 
even yet undiscovered mechanisms. 



lO 




10^ 10^ 10^ 

V, (km s-i) 



Fig. 2 Baryonic mass Mi, vs. circular velocity Vc in dwarf (squares) and spiral (circles) galaxies 
1891 and in clusters of galaxies (triangles; 11191 ). The blue dashed line is the fit to the spiral galaxies 
alone Mt = 50(V(./kms^')*M,.,. The red dotted line is the simplest standard model expectation if 
all the baryons in each DM halo are identified. Courtesy of S. McGaugh. Reproduced from 1901 . 
with permission. 

Tully and Fisher 01451 showed that, in disk galaxies, the maximum rotation ve- 
locity Vrot is proportional to the galaxy luminosity. MOND predicts that at large radii 

/oq ~ flAf = GMjj/r^ (equation[TJ and from a = v-^^Jrwe obtain 

vt^^GMhao. (3) 
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where is the baryonic mass. Figure|2]shows the measured baryonic mass M/, as a 
function of the measured circular velocity Vc for astrophysical systems on different 
scales, from dwarf galaxies (left bottom corner) to clusters of galaxies (upper right 
corner). Vc is a characteristic circular velocity measured in the outer region of the 
system 1901 . The blue dashed line shows the relation Mt, = 50(Vc/km s^ ' . The 
MOND prediction clearly agrees with both the observed slope and normalization. In 
addition, the observed spread is consistent with the uncertainties. This result implies 
that the relation v^o, = GMyao holds exactly. 

In the standard model, disk galaxies are embedded in DM halos whose aver- 
age density within their virial radius is basically independent of the halo mass. It 
is thus usual to define r2oo as the radius within which the average density is 200 
times the critical density of the Universe. It follows that for M200, the mass within 

r2oo, we have M200 °^ '"200 ™^ "^^e circular velocity Vc = (GM2oo/'"20o)'''^ °^ '^Im- 
This relation is shown as the red dotted line in Figure |2] . The circular velocity is 
not necessarily identical to the disk rotation velocity Vrot, because of the complex 
interplay between the merger history of the DM halo and the star formation history 
and energy feedback of the galaxy. For example, we could easily recover the correct 
slope y^Qj 0= Mh, by assuming that luminosity traces the total mass (L °^ Mh °^ Mjot) 
and that the density and scale height of the galaxy disk is roughly constant in disk 
galaxies; this latter assumption implies Mioi °= where R is the disk size. From 
^rot ^ GMtot/R, we correctly derive v^^f °^ Mh- A rigorous comparison between ob- 
servations and the standard model is not trivial |,137J . but recent analyses, where 
properly balanced contributions of the various physical and observational effects are 
carefully blended, seem to bring the ACDM Tully-Fisher relation in better agree- 
ment with observations 1381 l42l . 

However, in the standard model, the predicted scatter remains larger than ob- 
served, because, unlike MOND, we do expect that the galaxy merger and star for- 
mation history mentioned above introduce an intrinsic scatter In addition, a priori, 
we do not have any reason to expect that (1) the observed relation extends over five 
orders of magnitude in baryonic mass, from dwarf galaxies to massive disk galaxies, 
again with basically no scatter, as shown in Figure |2] and that (2) the LSB galaxies, 
that should presumably have a star formation efficiency lower than normal galax- 
ies 1150111461 . also perfectly fit into this relation. This latter result was originally 
predicted by Milgrom ll95l . fifteen years before the measurements of the rotation 
curves of LSB galaxies l92l . 

MOND goes beyond the description of global properties of disk galaxies. Fig- 
ure[3]shows the observed rotation curves of a HSB galaxy and a LSB galaxy (open 
symbols with error bars). The rotation curves expected in Newtonian gravity from 
the distribution of baryonic matter (black solid lines) severely underestimate the ob- 
servations. The underestimate increases with distance R from the galaxy center and 
is larger for the LSB galaxy. In the standard model, this observation is explained 
by increasing the DM contribution with increasing R and decreasing galaxy lumi- 
nosity. This solution is usually motivated by a lower star formation efficiencies at 
larger radii, as suggested by extensive surveys of neutral hydrogen in nearby galax- 
ies l8ni2T1. 
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Fig. 3 Examples of MOND rotation curve fits of a HSB galaxy (NGC 6946, top panel) and a LSB 
galaxy (NGC 1560, bottom panel). The black lines show the Newtonian rotation curve expected 
from the observed distribution of stars and gas. The blue lines are the MOND fits with best fit 
stellar mass-to-light ratios in the /sT-band 0.37 Mq/L,, (NGC 6946) and 0. 18 Mj/L v: (NGC 1560). 
Courtesy of B. Famaey and S. McGaugh. Reproduced from 1481 . with permission. 



What is more remarkable is that the small scale features of the rotation curves 
mirror the distribution of the baryonic matter in the disk. This characteristic is shared 
by the rotation curves expected in Newtonian gravity from the distribution of bary- 
onic matter (black soUd lines in Figure |3}- In the standard model, this property of 
the total rotation curve is unexpected, because this rotation curve should be mostly 
determined by the dynamically dominant DM distribution and the baryonic distribu- 
tion should play a very minor role. On the contrary, MOND describes the observa- 
tions with impressive accuracy (blue solid lines), including the small scale features 
of the curves. 

The MOND curves are derived with a single free parameter, the stellar mass-to- 
light ratio of the disk, in addition to the distance to the galaxy required to determine 
its length and hence acceleration scales; the best fit values of the mass-to-light ra- 
tios are in perfect agreement with values derived from stellar population synthesis 
models. Moreover, redder galaxies require larger mass-to- light ratios, as expected on 
completely independent astrophysical grounds lfT6l . This limited freedom of MOND 
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should be compared with the standard model that requires two parameters for the 
DM halo of each galaxy and a global mass-to-light ratio of the galaxy that depends 
on R and is unrelated to the mass-to-light ratios of the stellar populations. 

The agreement shown in Figure |3]is common to 75 nearby galaxies to date 14811 . 
Among these, an interesting case is NGC 7814: this galaxy is almost perfectly edge- 
on and the uncertainties on the rotation curve deriving from the disk inclination are 
negligible; in addition, the distance to the galaxy is accurate to 5% and is basically 
not a free parameter any longer This galaxy provides a stringent test for MOND: 
high-quality infrared photometric observations [!54| enable, for the first time, the 
construction, from an accurate bulge-disk decomposition, of a three-dimensional 
model of the galaxy whose gravitational potential is inferred by numerically solving 
the MONDian Poisson equation. The comparison of the model rotation curve with 
the observed one allows the derivation of the mass-to-light ratios for both the disk 
and the bulge components. Both ratios are found to be in excellent agreement with 
the expected values [lOj. 

The data expected in the upcoming GAIA mission and other future surveys will 
provide unprecedented possibilities to test MOND with the Milky Way dynamics 
ll48l . Pioneering work with current data has already shown that the rotation curve 
and the surface density of the inner disk of the Milky Way are fully consistent with 
each other within the MOND framework ||9T1 1471 . Similarly, the escape velocity 
from the solar neighbourhood agrees with current estimates if the external gravita- 
tional field in which the Milky Way is embedded is a ~ lO^^ao: this value is indeed 
compatible with the actual field 1461 . An additional test involves the velocity ellip- 
soid tilt angle within the meridional galactic plane. The angles expected in MOND 
and in the standard Newtonian gravity with DM agree with each other and with 
observations at galactic heights z = 1 kpc; however, the discrepancy between the 
predicted angles in the two models increases with z and the measure of the velocity 
ellipsoid tilt angle will thus be a relevant test to discriminate between the models 
iQl. 



3.2 Elliptical Galaxies 

The role played by the acceleration scale is also evident in elliptical galaxies, dwarf 
spheroidal galaxies (dSphs) and globular clusters. 

In elliptical galaxies, that are pressure supported systems, Faber and Jackson l45l 
observed a relation similar to the Tully-Fisher relation that is valid for disk galax- 
ies: the galaxy luminosity L correlates with the stellar velocity dispersion a in the 
galaxy's central region, according to the power law L °^ (7^. If we assume that ellip- 
ticals are isothermal spheres that, in MOND, have finite mass with asymptotically 
decreasing density p °^ r", with a ~ d\np/d\nr = —4, we find a'^ = GMi,ao/a^ 
l97ll . Unlike the Tully-Fisher relation, the observed Faber-Jackson relation is not 
expected to be exact in MOND, because ellipticals are not strictly isothermal and 
their velocity field is not isotropic; the velocity anisotropy parameter must actually 
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vary with radius to match the observed Fundamental Plane of ellipticals 112411 . The 
MOND Fundamental Plane is actually slightly tilted compared to observations, but 
this problem might be removed by including the external field effect 1281 . 

Elliptical galaxies pose some challenges to the standard model, because they 
should be embedded, like disk galaxies, in DM halos. X-ray emitting hot gas coro- 
nae are expected signatures of a DM halo and are indeed observed in many early 
type galaxies |ll3j, including isolated ellipticals, like NGC1521 l65i . However, 
there are cases that are unexpected in the DM paradigm: accurate observations, 
based on planetary nebulae, of the kinematics of the outer parts of three ordinary 
ellipticals show very little evidence, if any, of the presence of DM 1121], but are in 
good agreement with MOND, because the large masses implied by the high surface 
brightnesses indicate that the gravitational field is in the Newtonian regime a > oq 
005111421 . 

The kinematics of the outskirts of ellipticals, where a < gq, can be probed with 
spectroscopic observations of their globular cluster systems: for example the galaxy 
NGC4636 in the Virgo cluster is surrounded by 460 globular clusters with measured 
velocity and this sample represents one of the largest currently available. The MON- 
Dian predictions agree with the kinematic data and NGC4636 also appears to fall 
onto the baryonic Tully-Fisher relation shown in Figure |2l [|128l| . 

The X-ray data available for NGC1521 and NGC720 offer a unique test of 
MOND in elliptical galaxies: these galaxies are not embedded in groups or clusters 
and the X-ray data extend to large radii. Similarly to disk galaxies, one can thus test 
MOND in the outskirts of galaxies where the gravitational acceleration due to the 
luminous matter is smaller than oq and the external field effect is negligible. MOND 
describes the distribution of the baryonic mass in these galaxies with mass-to-Ught 
ratios fully consistent with stellar population synthesis models 01041 . 



3.3 Dwarf Spheroidals 

At the low mass end of the galaxy mass function, dSphs also pose a challenge to the 
standard model (see 1721 for a recent review of dSphs in MOND). dSphs have low 
surface brightnesses and, according to MOND, based on what is shown in Figure [T] 
for disk galaxies, Milgrom predicted that they should be DM dominated and have 
mass discrepancies larger than ten when analyzed with Newtonian gravity 1,95] . This 
prediction was impressively confirmed when the first measures of the stellar velocity 
dispersion in the central regions of these galaxies were available a decade later |98 1. 

More recently, intense observational programs provided velocity dispersion pro- 
files of the dSphs orbiting the Milky Way II 1491 . This piece of information enables 
the estimate of the combination of the mass profiles of the galaxies with the pro- 
files of their velocity anisotropy parameter The most recent detailed dynamical 
analysis II129I confirms that the mass-to-light ratios in the V band are in the range 
1—3 Mq/Lq, and are therefore consistent with stellar population synthesis models. 
This analysis solved an open issue raised earlier [2]: unbound stars can contaminate 
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the velocity dispersion and artificially inflate the estimate of the mass-to-light ratio. 
When the unbound stars are properly removed with the caustic technique ||40lll30L 
Sculptor and Sextans do indeed show the sensible mass-to-light ratios 1.8 Mq/Lq 
and 2.7 Mq/Lq respectively, whereas Carina still shows a too large ~ 6 MqJLq. 

This discrepancy might originate from (1) the uncertainties of the luminosity dis- 
tribution that is challenging to estimate accurately enough in LSB galaxies, and (2) 
the ellipticity of Carina that significantly departs from the spherical symmetry as- 
sumed to derive the mass-to-light ratio. However, this discrepant mass-to-light ratio 
might have a deeper origin due to the specific feature of MOND that we mentioned 
in Sect. |2l the external gravitational field effect: a star in a dSph, that is a satellite 
of the Milky Way, moves according to both the dSph mass and the gravitational 
field exerted by the Milky Way; only if this latter external acceleration is negli- 
gible compared to the acceleration internal to the dSph is the dSph mass derived 
from the stellar velocity dispersion accurate. Carina is one of the least luminous, 
and presumably least massive, dSphs and one of the closest to the Milky Way. It is 
therefore reasonable to suspect that the Milky Way's gravitational field can play a 
role in inflating the velocity dispersion of this dSph. Although this suggestion still 
awaits a quantitative confirmation, this same effect appears to be responsible for the 
deviation of some dSphs from the expected Tully-Fisher relation shown in Figure|2] 

An additional piece of evidence, which is problematic for the standard model, 
is the phase-space distribution of the dSphs that are satellites of the Milky Way. 
These dSphs are distributed over an extended thin disk whose thickness is between 
10 and 30 kpc and radius ~ 200 kpc. An invoked solution is that these dSphs fell 
into the MW halo as a small group of galaxies who kept their orbits correlated ll82l . 
However, recent measures of the dSph proper motions indicate that this scenario is 
untenable because, according to these measures, four dSphs must have fallen in at 
least 5 Gyr ago lfT2l and A^-body simulations of the standard model show that the 
orbit correlation cannot be preserved for such a long time [|70||| . 

In MOND, dSphs may form as tidal debris during close encounters of large galax- 
ies ll76l and the orbit correlation would thus be a natural consequence of this for- 
mation process 111 14111 15l . The formation of tidal dwarf galaxies has been observed 
in interacting galaxies for the last twenty years [41], since the first detection in the 
Antennae |106|. The observed stellar velocities of these systems agree with MOND 
in the majority of the systems using sensible values of the two available free, but 
constrained, parameters, the mass-to-light ratio and the velocity anisotropy param- 
eter /3 moil l56l . In the standard model these systems are somewhat challenging 
because, in this case, the observed velocities require a factor of a few more mass 
than the observed luminous mass, despite the fact that there is no physical reason 
for these low surface brightness tidal dwarfs to drag large amounts of DM 1221 . 

The formation of dwarf galaxies as tidal debris is also likely to solve the missing 
satellite problem of the standard model that predicts a factor of ten more satellites 
in the halos of large galaxies than is actually observed 010711 . In MOND, the rate of 
galaxy encounters is small enough that it might provide the right number of satel- 
lites with the correct dynamics. On the contrary, in the standard model, a conspiracy 
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of processes regulating the star formation efficiency is required so that most DM 
satelUtes form no stars ||62) . In addition, the possible supersonic relative velocity 
between baryons and DM before reionization might be responsible for the inhibi- 
tion of the formation of half the expected luminous satellites fTSl. If we take into 
account the ultra-faint galaxies surrounding the Milky Way, the problem can also be 
partly alleviated, but there would still be a factor of four too few dwarfs 11341 . In 
addition these ultra-faint galaxies show indications of tidal disruption, although one 
would expect that, with mass-to-light ratios of the order of 1000 Mq /Lq, the large 
DM halo within which they are embedded should be sufficient to screen the stellar 
components from the external tidal field of the Milky Way. The properties of dwarf 
galaxies thus remain challenging for the DM paradigm ||25] l64l l75l . 



3.4 Globular Clusters 

The good agreement between the expected dSph dynamics in MOND and observa- 
tions is more intriguing when we consider globular clusters. These stellar systems 
and dSphs roughly have the same baryonic mass but different surface brightnesses, 
stellar populations and ages. In the standard model, globular clusters are devoid of 
DM, whereas dSphs are the cosmic structure with the largest fraction of DM, with 
mass-to-light ratios in the range ~ 10- 1000 Mq/Lq II 1491 11481 . These discrep- 
ancies between globular clusters and dSphs are solved by invoking low formation 
efficiencies in low mass DM halos 1741 l73l and two completely different formation 
processes for the two kinds of systems. On the contrary, in MOND, the different 
observed internal velocities is exactly what we expect from the different surface 
brightnesses (Figure [T) [55|. 

Nevertheless there is a system that requires careful consideration. NGC 2419 is a 
globular cluster in the outer halo of the Milky Way with low enough surface bright- 
ness to be in the MOND regime. It is far enough from the gravitational field of the 
Galaxy that the MOND external field effect might not play a relevant role. Similarly 
to the analyses of the dSphs, it is possible to estimate the mass-to-light ratio of the 
cluster from measurements of the projected velocity dispersion of the stars ll66l . In 
MOND, the mass-to-light ratio required to describe the stellar kinematics is a factor 
of a few lower than expected (the opposite case as for the Carina dSph). Invoking 
a variable velocity anisotropy does not help improve the comparison between the 
MOND fit and the data, but non-isothermal polytropic models seem to provide a 
MONDian description of the kinematic and photometric observations of this cluster 
II126II1271 . An additional solution is that the life time of the globular cluster is long 
enough that mass segregation has already taken place. In this case, the distribution 
of RGB and luminous upper main sequence stars used in the Jeans analysis is ex- 
pected to be more centrally concentrated than the other less massive stars of the 
clusters; therefore the true gravitational potential of the cluster might be different 
from the gravitational potential actually derived with the Jeans equation. Another 
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possibility is that globulars have non-standard initial stellar mass functions. More 
data and more testable globular clusters are clearly called for 

Another system that has seemed to be challenging for MOND is the five globular 
clusters orbiting the Fornax dSph. Fornax is the most luminous dwarf spheroidal (by 
about a factor of ten) and is the only classical dwarf of the Milky Way with a system 
of globulars. It was suggested that the surprisingly stronger dynamical friction due 
to the dwarf galaxy's low density component of stars in MOND, relative to the 
dynamical friction from the more dense DM halo in Newtonian gravity 041 . would 
cause the globular clusters to lose their orbital angular momentum in a period of time 
much shorter than the Hubble time and create a stellar nucleus in Fornax which is 
not observed 11221 l36l . However, Angus and Diaferio [[41 used an orbit integrator 
with accurate mass models of Fornax both in MOND and Newtonian dynamics to 
convincingly demonstrate that the situation is relatively easy to explain: the globular 
clusters can orbit for a Hubble time as long as their orbits start near the tidal radius. 
This solution does not apply to the most massive globular cluster; however, this 
globular cluster is a statistic of one and could have a sizeable line of sight distance 
where the dynamical friction is negligible. 



3.5 Groups and Clusters of Galaxies 

The impressive agreement between observations and MOND predictions on the 
scale of galaxies, based on the introduction of the acceleration scale oq, is not shared 
by the galaxy cluster data Ill39ll57l . Clusters do not perfectly obey the relation be- 
tween circular velocity and baryonic mass of equation (|3) (green triangles in Figure 
|2]i but seem to require more mass than actually observed. In fact, in the core of 
clusters a > qq and the luminous matter should be enough to describe the observed 
dynamics. The evidence suggests the opposite is true, because the amount of ob- 
served mass is a factor of two too small. The shortage of mass, which is confined to 
the central regions of clusters and progressively disappears in the cluster outskirts 
Q, is not as large as the factor of five or larger we have in the standard model, but 
it clearly poses a challenge to MOND. 

An equivalent, more conventional way to look at Figure |2]for clusters is to con- 
sider the relation between the cluster mass and the X-ray temperature of the hot in- 
tracluster gas. According to the relation between mass and circular velocity, MOND 
makes a clear prediction for the mass-temperature relation. By assuming the same 
argument used for the Faber-Jackson relation in elliptical galaxies, the relation be- 
tween the velocity dispersion CT of the galaxies in a cluster and the baryonic mass M/, 
of the cluster is = GM^ao / 0!^, where a is the logarithmic slope of the mass den- 
sity profile. The temperature T of the intracluster medium is a measure of the kinetic 
energy of the galaxies and T °^ . Therefore, in MOND, we have M\, <x T^. Figure 
|4]shows that the data agree with this scaling relation but not with the normalization, 
because of the mass discrepancy in the cluster cores mentioned above. 
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Fig. 4 The baryonic mass-X-ray temperature relation for rich clusters (gray tiiangles 11251 ) and 
groups of galaxies (green triangles (?]). The solid line is the MOND prediction Mb <x T^. Courtesy 
of B. Famaey and S. McGaugh. Reproduced from |48:|, with permission. 



This relation is not at odds with the standard model relation M oc r^/^H which is 
known to agree with observations [il47( . for two reasons. First, the M °^ T^l^ relation 
is between the X-ray temperature and the total cluster mass, which includes DM, 
rather than the baryonic mass alone that we have in the MOND Mj, °^ relation. 
Second, the mass M is inferred by assuming valid Newtonian gravity whereas Mj, 
shown in Figure|4]is derived with MOND \T\. Nevertheless, we can neglect the latter 
reason and still see that the two scaling relations are roughly consistent with each 
other In fact, in the standard model, we can write M = Mi, / fb°^ T^l'^, where //, is the 
baryon fraction of the cluster total mass. Therefore, the observed MONDian M/, 

should imply fb^T^^^- This expectation is broadly consistent with observations: 
by combining X-ray Chandra groups and HIFLUGCS 111511 over the temperature 
range [0.6, 15] keV, Eckmiller et al. Ell find ~ r0.79±o.09 ^^^^.^ |,y Umiting the 
sample to clusters with T in the range [1,15] keV, the slope is 0.83 ±0.42, but it 
appears to be shallower at smaller radii. 



' In the standard model, the cluster virial mass scales as M^-„ x p^.{z)Ac{z)R^ , where R is the cluster 
size in proper units (not comoving), (z) — (z) / (8;rG) is the critical density of the Universe, 
and Ac (z) is the cluster density in units of (z) . A widely used approximation is 



where w = n„,{z) - 1 |27|. Now, p,.{z) scales with redshift z as p,.{z) E^{z) = 0,„{l + zf + 
(1 - n,„ - i2A)(l -|-i2A. The cluster size thus scales as R mI[^Ac^^^{z)E-^/^{z), and the 
temperature as T <x My^/R <x mV^^'''^ {z)E^/^ (z). 
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The apparent mass discrepancy in X-ray bright groups and clusters in MOND Q, 
clearly shown by the difference between the observed and predicted normalizations 
of the mass-temperature relation shown in Figure 2] does not seem to appear in 
groups of galaxies where the X-ray emission is negligible or absent 199111001 . This 
piece of evidence might suggest that the mass discrepancy in MOND can be solved 
by assuming the presence of undetected baryonic mass in the form of cold, dense 
gas clouds that can form and be stable because of the presence of the ionised hot gas 

IMI. 

A different solution is that we might well have some DM that clusters on cluster 
scales but not on galaxy scale. The piece of evidence that is usually claimed to be 
difficult to interpret without assuming the existence of some form of DM comes 
from colliding galaxy clusters. During the collision of two clusters, the galaxies 
and the two halos of DM, if DM exists and is collisionless, keep their momentum, 
whereas the intracluster medium is shock heated and slows down. The DM and 
galaxy components thus separate from the gas component. This separation is clearly 
visibile by combining X-ray and optical images of individual colliding clusters: 
for example, the so-called Bullet cluster IE 0657-558 135], and the cluster MACS 
J0024.4-1222 [|26l . 

In the absence of DM, most of the matter resides in the gas rather than in the 
galaxies. The only way to measure how the total mass is distributed, is to derive a 
map of the gravitational potential with weak gravitational lensing. In GR, this pro- 
cedure is now standard and in the two systems mentioned above, the mass appears 
to be concentrated where the galaxies are and not where the gas is: this is in striking 
agreement with the presence of some form of collisionless DM dominating the mass 
content of the cluster. 

Unfortunately, a similar gravitational lensing analysis is not trivial in alternative 
theories of gravity where the effect of gravity on the light path is not treated prop- 
erly. MOND belongs to this group of theories because it is a classical theory and 
gravitational lensing can only be described by resorting to one of the different co- 
variant extensions of MOND. We will describe the results of this approach in Sect. 
|4] Here, we wish to emphasize that the location of the gravitational peaks is not an 
observational fact, but derives from the assumption of GR's validity. A priori it is 
possible that the non-linearity of a gravity theory alternative to GR is sufficient to 
mimic the observed gravitational lensing distortion with a mass distribution differ- 
ent from that required by GR. For example, in the covariant MOND theory named 
TeVeS, particular matter distributions can yield non-zero convergence on sky po- 
sitions where there is no projected mass ||9]. Weyl gravity provides another clear 
non-MONDian example of this phenomenology 11 16111 ITl . 

There is a final cautionary comment that cannot be omitted: the interpretation of 
the observed properties of colliding clusters is far from being clear, because these 
systems can also be challenging for the standard model. The Abell cluster A520 has 
the same morphology of the Bullet cluster and MACS J0024.4-1222: two clouds of 
galaxies on the opposite sides of a cloud of hot gas 18411671 . However, the GR weak 
lensing analysis indicates that, in addition to the two gravitational potential peaks 
at the location of the galaxies, there is a significant peak where the gas is. This 
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peak is difficult to explain in the standard model, because it requires the existence 
of a massive DM halo devoid of galaxies, but associated with the cluster hot gas. A 
number of possible solutions have been suggested, but none of them appears to be 
fully convincing 1671 . 

In addition, colliding clusters appear to have relative velocities that are unlikely 
in the standard model: a system like the Bullet cluster, with a relative velocity de- 
rived from the Shockwave of ^ 4700 km s^' ||35l . requires an initial infall velocity of 
the two clusters of ^ 3000 km s^ ' 1851 that has a probability smaller than 3 .6 x 10^^ 
to occur in ACDM 1801 . On the contrary, the enhanced intensity of MOND grav- 
ity may naturally produce these large relative velocities [6 5|. Similarly, the co- 
herent motion of galaxy clusters on large scales, measured with a technique based 
on the kinematic Sunyaev-Zeldovich effect produced by Compton scattering of the 
CMB photons, appears to be challenging for the standard gravitational instability 
paradigm: these bulk flows are a factor of five larger than predicted by the stan- 
dard model and might require a modification of the theory of gravity, among other 
possible solutions 1681 . 



4 TeVeS and Gravitational Lensing 

MOND, as described in Sect. |2l is a classical empirical law that we would like to 
recover as the weak field limit of a covariant theory. This theory should contain GR, 
that excellently describes the gravitational phenomenology of the solar system. A 
covariant theory for MOND is required if we wish to build a cosmological model 
and describe the gravitational lensing phenomenology within the MONDian frame- 
work. A covariant theory is thus essential for the validation of MOND. 

The attempts to build a covariant theory containing MOND are numerous. A 
recent overview of these attempts is provided by Famaey and McGaugh |48J. Here, 
we briefly outline one of them, namely TeVeS (see lT4l for a recent review), whose 
Lagrangian contains a time-like vector field and a scalar field in addition to the 
tensor field representing the metric, hence the name Te(nsor)Ve(ctor)S(calar). Halle 
et al. l63l have shown that a very general Lagrangian for a decaying vector field 
unifies most of the popular gravity theories, including quintessence, f{R), Einstein- 
Aether theories, and TeVeS itself, among others. We stress that TeVeS by no means 
is the only possible covariant theory that yields MOND in the weak field limit. 
Moreover, it is clear that any problem or failure belonging to one of these covariant 
theories are not necessarily problems or failures of MOND. 

TeVeS was first proposed by Bekenstein ||13| . At the classical level, MOND in- 
troduces an acceleration scale. This issue poses an immediate problem for building 
a generally covariant theory, because the acceleration scale is played by the affine 
connection f^y, that involves the first derivatives of the metric tensor gfjy, and 
is not a tensor As we mentioned in Sect. |2] one way to bypass this problem is to 
distinguish between the Einstein metric g/jy, that enters the Einstein-Hilbert action, 
and the geodesic metric gfjy, that enters the matter action. In GR, the two metrics 
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coincide. In TeVeS, the two metrics are related by the equatior0 II123I 

lMv=e-2%v-2sinh(2(p)C/^f/v, (5) 

where ^ is a scalar field and f/y^ is a normalized vector field with g^^U^iUy = — 1. 
Both fields are dynamical. Therefore, the TeVeS action is the sum of three terms: 
the standard Einstein-Hilbert action, the action term for the scalar field 

= "2^ / ^i''^^h"^M.Ii)i~8)"^^'^ ' (6) 

where ^ is an arbitrary positive function, A: is a dimensionless coupling constant, £ 
a constant scale length, and /i"^ = g"-^ — g"^g^^Ufj,Uv, and the action term for the 
vector field 

% = - ^ /[^^"''^'''t/[a,M]^/3,vl +^(^"%a;/3 -2A fe^'Xf^v + 1 )] (-^) ^/'d 

(7) 

where the square brackets indicate antisymmetrisation, K and K are dimensionless 
coupling costants and A is a Lagrange multiplier to guarantee the normalization of 

TeVeS violates the local Lorentz invariance, because at each point in spacetime 
there is a preferred frame in which the time coordinate aligns with U^. The vio- 
lation of Lorentz invariance derives from the invalidity of the Strong Equivalence 
Principle anticipated in Sect. |2] Clearly this violation has to be smaller than current 
experimental bounds. 

In the limit K,K,l/£ ^ 0, with k ^ £^^/^, for quasi-static systems and homo- 
geneous cosmology, TeVeS corresponds to GR. To recover MOND in the non- 
relativistic ultra- weak field limit (equation [T}, we need to choose the function JF. 
The arbitrariness of the function ^ makes TeVeS a family of models, rather than 
a single model. The function ^ proposed by Bekenstein ifTSl yields the MOND 
acceleration scale 



4. 



As we have just mentioned, we recover GR in the limit K,K ^0, whereas to recover 
Newtonian gravity in the weak field limit £ — > oo suffices. Therefore, in principle, 
Te VeS and GR might differ in the strong field regime [79 1 . 

In the weak field limit, and quasi-static system, the geodesic metric becomes 

l^vdjc^'djc*' = -(1 +2<?)df2 + (1 +2f )5,7djcW . (9) 

This metric is formally identical to GR, where <P ^ ~ (j)^ and (pN is the New- 
tonian gravitational potential. In this limit, in TeVeS we also have = — f , but 
~ E^N + (p, where E ~ {1 — K /2)^^e^^f'-' ^ 1, with (pc the asymptotic boundary 



* In this section, we use units where the speed of light c = 1 . 
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value of (p. It appears clear that the scalar field (p plays the role that DM plays in the 
standard model. 

Equation (|9]l enables us to derive the gravitational lensing equations in TeVeS. In 
Sect. |3]we have seen that MOND excellently describes the dynamics of galaxies but 
requires some additional DM on the scales of clusters. One can thus anticipate that 
the description of gravitational lensing with TeVeS fits the observations on galaxy 
scales without the need of any matter in addition to the observed baryonic matter, 
but it does not do as well on cluster scales. 

In fact, the TeVeS gravitational lensing equations applied to the strong lens- 
ing regime on galaxy scales reproduce the lensing morphology with the observed 
baryonic matter alone, both in simple spherically symmetric models of the lens 
II32I I152II and in models departing from spherical symmetry I1131II . These analy- 
ses self-consistently include TeVeS cosmology, because of the cosmological dis- 
tances of the lenses: it is easy to get erroneous results if one uses hybrid theories 
hke MOND combined with GR rather than TeVeS itself [531. In addition, claimed 
inadequacies of TeVeS lf86l may be easily healed by assuming different, but still per- 
fectly reasonable, mass models for the lens ||33l . When analysed properly, to date no 
lensing system on galaxy scale appears to be problematic, including lensed quasars 
£31 , 30 33j. TeVeS even returns a measure of the Hubble constant consistent with 
other independent estimates when its time-delay formula is applied to lensed vari- 
able quasars 01411 . Nevertheless, the debate on the adequacy of TeVeS to describe 
lensing data on galaxy scales is still open and lively ll52l . 

Lensing of clusters returns the mass discrepancy that one encounters with the 
analysis of galaxy kinematics or X-ray emission. Because of the complexity of the 
covariant theories of MOND, the lensing analysis usually assumes spherical symme- 
try and quasi-stationarity. Relatively relaxed clusters, where these assumptions are 
substantially valid, show the expected mass discrepancies 11381110811 . DM halos are 
required, and it was suggested that DM in MOND could be made of neutrinos 01251 . 
However, the neutrino phase-space density must be smaller than the Tremaine-Gunn 
limit we describe below in Sect. 15.11 The combination of strong and weak gravita- 
tional lensing data shows that the DM central density of MONDian clusters is larger 
than the Tremaine-Gunn limit for neutrino masses smaller than 7 eV 11081 : this 
mass exceeds by more than a factor of three the experimental upper limits on the 
neutrino mass and this result suggests that neutrinos are inadequate MONDian DM 
particles. 

The Bullet cluster is usually held as the definitive proof of the existence of colli- 
sionless DM ll35l . Unfortunately, it is a difficult lens to model in MOND covariant 
theories because of its large departure from the assumptions mentioned above of 
spherical symmetry and quasi-stationarity. Figure |5]shows the results of one of the 
first TeVeS analyses of the Bullet cluster ifTTI : the convergence map is computed by 
assuming four spherically symmetric mass distributions at the location of the two 
galaxy distributions and the two X-ray emitting gas clouds. The model agrees with 
observations if the galaxies are embedded in additional halos of collisionless mat- 
ter Angus et al. jTT] concluded that massive neutrinos with 2 eV mass accounting 
for two thirds of the total mass of the system is sufficient to reproduce the lens- 
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Fig. 5 The solid black contours show the MOND-ian convergence map of the Bullet clus- 
ter nil . The dotted red contours show the GR convergence map |35| . The contour levels are 
[0.37,0.30,0.23,0.16]. The red stars indicate the centres of the four potentials used. The blue 
dashed lines show the contours of surface density [4.8, 7.2] x 10" M;., pc^- for the MOND standard 
fl function. In the green shaded region the matter density is larger than 1.8 x 10^"* Mq pc^-' and 
indicate the clustering of 2 eV neutrinos. Mset: The surface density of the gas in the Bullet clus- 
ter predicted by a collisionless matter subtraction method for the standard /i-function described 
in [ni. The contour levels are [30,50,80, 100,200,300] M,,,pc--. The origin is [06'' 58"' 24'. 3 8, 
-55°56'.321. Reproduced from fTTI. 



ing signal. This amount of additional mass agrees with the results of the dynamical 
analysis of other clusters in MOND reviewed in Sect. 13.51 but, it became clear later 
iTTl llOSI . as we mention earlier, that the conclusion about the 2 eV mass neutrinos is 
erroneous. The non-linearity of TeVeS does not seem to be conducive to removing 
the demand for DM in the Bullet cluster ifSOl . On the other hand, the effect of the 
external gravitational field, that plays a role in the internal dynamics of dSphs, or 
the non-trivial features that the additional TeVeS fields can induce in the lensing 
phenomena remain to be investigated ISTl . 



5 MOND and Sterile Neutrinos 

MOND assumes the existence of an acceleration scale below which an ultra-weak 
field limit of the theory of gravity sets in. From this very simple ansatz, MOND 
describes with impressive success the dynamics of galaxies and smaller systems. 
On the scale of galaxy clusters, MOND partly fails (Figure |2) and the observed 
kinematics require either a new gravity law or the existence of some baryonic or 
non-baryonic DM. This failure is bound to be shared by any covariant theory, like 
TeVeS, that is conceived to yield MOND in the ultra-weak field limit. 
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For most cosmologists, this failure implies the end of MOND as a successful 
description of the Universe. However, the DM paradigm currently is far from being 
satisfactory on galactic scales JTS). Therefore, if we are interested in finding a the- 
ory that describes the Universe with the minimum number of assumptions, we have 
to consider the possibility that MOND can indeed be a valid description of the ob- 
served phenomenology and look for possible solutions of its shortcomings on larger 
scales. 

Similarly to the standard model, a natural solution is the introduction of some 
form of DM. However, we might have an advantage over the standard model. In a 
MONDian model with DM, DM has to be hot enough to freely stream out of galactic 
systems, to preserve the excellent description of the galactic dynamics without DM, 
but cold enough to cluster on the scale of galaxy clusters. Unlike the hypothetical 
cold DM particles, we know that an elementary particle that can play the role of hot 
DM does exist: the neutrino. 

The properties of neutrinos are currently constrained by various experimental re- 
sults. In 1979, Tremaine and Gunn 01431 . by considering the maximum mass density 
that DM halos made of light leptons can reach, set a lower limit of ~' 1 MeV to the 
mass of the light leptons that can make the DM halos of galaxies and clusters of 
galaxies. Below, we briefly review this argument for its relevance to the subsequent 
discussion. 



5.1 The Tremaine-Gunn Limit 

Neutrinos are collisionless particles, and, according to Liouville's theorem, the 
phase-space fluid they form is incompressible. In practice, if neutrinos make the DM 
in self-gravitating systems, like galaxy clusters, this theorem sets an upper limit to 
the observable coarse-grained phase space density. 

If self-gravitating systems form by violent relaxation, the neutrino coarse-grained 
phase-space density is /(x, v) =/(e) =/o{l +exp[j3(e — x)]}^^ where e = v^/2 + 
^(x), 0(x) is the gravitational potential, (J^(x) = 1//3 is the ID velocity dispersion, 
X is the neutrino chemical potential, and /o = gyniyh^^ is the mass phase-space 
density of an occupied microcell; gv = 2 is the number of degrees of freedom, which 
includes the anti-particles, niy is the neutrino mass, and h the Planck constant. 

For a non-degenerate neutrino fluid, we have /(e) <C fo which implies ^{^ — 
X)^^ and the phase-space density must be smaller than the Maxwell-Boltzmann 
distribution /mb(£) = /oexp[— v^/2(7^(x)]. Therefore, for the neutrino mass den- 
sity Pv(x), we must have Pv(x) < An fQ°° fuE{£)dv = fo\2nG^{x)f'l^ . This re- 
lation implies that max{/} = /o > Pv(x)[2;rc7^(x)]^^/^. However, clusters form 
from the relic neutrino background that has the Fermi distribution f{p) = /o[l + 
exp(/7c/A:Br)]^', with p the momentum and T = 1.95 K the neutrino tempera- 
ture today. The initial maximum coarse-grained phase-space density is therefore 
max{/} = f{p = 0) = /o/2. According to Liouville's theorem, this upper limit can- 
not increase, and we must thus have 
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p,{x)<^[27za\x)]y\ (10) 

For a fully degenerate gas, all the microcells with v < viim are occupied and 
/(e) = /o, whereas all the microcells with v > vnn, are empty and /(e) = 0. 
Therefore, in this case, the neutrino mass density is Pv(x) < fo4n /q''""v^c/v = 
/o47rv^jj^/3. By considering only bound states with e = v^/2 + 0(x) < 0, we have 

< 2\(j)\ = v^i^, and thus Pv(x) < /o4;r|2(/)p/^/3. By assuming 3(J^ = —0, we 
obtain Pv(x) < /o[27i:o'^(x)]^/^4(3/7r)'/^, which implies that the maximum phase- 
space density is max{/} = fo > pv(x)[27r(7^(x)]^^/^(3/7i:)^'/^/4. However, apply- 
ing Liouville's theorem again with the initial maximum phase-space density /o/2, 
we obtain the more stringent upper limit to the neutrino density 



Pv(x)<^[2;r(T2(x)]3/24f- ) , (11) 



1/2 

2 -yn^ 

which is 4(3/7r)'/^ « 3.91 times larger than the density upper limit obtained in the 
non-degenerate case 1781 11441 . Therefore the non-degenerate case yields the most 
restrictive upper limit and is usually called the Tremaine-Gunn limit. 
We can write this limit in astrophysical interesting units as 



Pv<2.16xl0^^-J [-) ^ iU) 

or, by considering the relation k^T / iinip = between temperature and velocity 
dispersion, with nip the proton mass and jj, = 0.6 the mean atomic weight for a fully 
ionized gas of solar abundance. 



5.2 The Role of an 11 eV Sterile Neutrino 

The mass contribution of the three ordinary (active) neutrinos of the standard model, 
Ve, Vfi, and Vt, could, in principle, solve the mass discrepancy in MONDian galaxy 
clusters. At the end of the nineties, laboratory experiments yielded an upper bound 
Umit to the mass of ordinary neutrinos of 2.2 eV lISTI . Therefore, neutrinos with 
mass 2 eV were proposed as DM in MONDian models 01251 . However, when 
their contributions to the properties of astrophysical sources, namely galaxy clus- 
ters, CMB and massive galaxies, are analyzed in details, the ordinary neutrinos with 
mass smaller than this limit are shown to be inadequate. 

Sanders 01251 and Pointecouteau and SiUc jl 181 pointed out the relevance of 
the Tremaine-Gunn limit in galaxy clusters in MOND. An extensive analysis of 
26 X-ray emitting groups and clusters of galaxies |j7] considered neutrinos with 



22 



Antonaldo Diaferio and Garry W. Angus 



the maximum mass allowed by the current upper hmits. These neutrinos have a 
Tremaine-Gunn limit that is at least a factor of two smaller than the DM density 
within 100 kpc that is required to describe the thermal properties of the intra-cluster 
medium in MOND. Therefore, either neutrinos are more massive, that is excluded 
by laboratory measurements, or they are not the major contribution to the MOND 
DM in clusters. 

The three species of ordinary neutrinos with mass in the range 1 — 2 eV also 
are problematic when we attempt to reproduce the CMB power spectrum. These 
neutrinos suppress the third peak by ~ 25% when compared to observations, be- 
cause of their free-streaming imposed by the Tremaine-Gunn limit ll88l . We can 
obtain a CMB power spectrum consistent with observations with these neutrinos 
in TeVeS if we substantially increase the MOND acceleration scale aq at the time 
of recombination II133I . However, this redshift dependence of oq remains unproved 
observationally. 

Finally, the amplitude of the weak lensing signal around luminous galaxies (L > 
IO^'L©), extracted from the Red-Sequence Cluster Survey and the Sloan Digital 
Sky Survey, suggests mass-to-light ratios larger than ^ 10 Mq/Lq in the MOND 
framework II140I . It thus follows that DM is also required on the scale of massive 
galaxies. This DM cannot be made of ordinary neutrinos that are not massive enough 
to cluster on these small scales. 

The inadequacies of ordinary neutrinos with any mass smaller than the current 
upper limits to properly describe this astrophysical phenomenology forces us to 
conclude that ordinary neutrinos as the MONDian DM are ruled out. 

MOND can still be viable if we resort to hot DM made of sterile neutrinos. Ster- 
ile neutrinos do not have standard weak interactions and are right-handed, unlike the 
three ordinary neutrinos, and are motivated by a number of anomalies observed in 
neutrino experiments (see JT] for an extensive review). For example, the existence 
of one sterile neutrino in addition to the three ordinary neutrinos can elegantly ex- 
plain the disappearence of electron neutrinos from the low energy beam measured 
in short-baseline neutrino oscillation experiments ll59] 16011771 1 10911 . 

If sterile neutrinos are more massive than ordinary neutrinos, they can have larger 
Tremaine-Gunn limits, and thus eventually solve the problems we mentioned above. 
The analysis of the MiniBooNE experiment data does indeed favour a sterile neu- 
trino mass in the 4—18 eV range liSSl . This mass range is inconsistent with cos- 
mological constraints in the ACDM framework that would require the mass to be 
smaller than 1 eV (see, e.g., Ii60l and references therein). However, this mass limit 
derives from the enhanced gravity on galactic scales that is required by the standard 
model but not by MOND. Therefore, in MOND, the 4 — 18 eV mass range might be 
perfectly reasonable. 

It is actually quite impressive that a universe with baryonic matter and one mas- 
sive sterile neutrino alone reproduces very well the observed CMB power spectrum 
if the sterile neutrino mass is niv, = 11 eV lO, a value that is fully consistent with 
the mass range inferred from the MiniBooNE data. In addition, the sterile neutrino 
mass must be within ~ 10% of the 1 1 eV value if we wish to keep the good match 
with the observed CMB power spectrum. This strong constraint is mostly due to the 
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fact that, in this case, the contribution of the sterile neutrinos to the density of the 
Universe is = 0.1 17. This value is comparable to the contribution of CDM in 

the standard model that we know to describe the CMB power spectrum very well. It 
is important to remind that this analysis assumes that the sterile neutrinos are fully 
thermalised at the time of decoupling; more massive sterile neutrinos are possible 
but they must not be fully thermalised j24l. 

Figure |6] shows that the CMB power spectra of the standard A CDM model and 
of a model with baryonic matter and an 1 1 eV sterile neutrino are indistinguish- 
able. Note that at the time of recombination, the average gravitational field is strong 
enough that the Universe is not in the MOND regime, with ciq kept constant to the 
present day value; therefore, the CMB power spectrum can be estimated with the 
standard theory of gravity. It is only at later times that the existence of such a mas- 
sive neutrino would be problematic for the formation of structure in the standard 
model. 




Fig. 6 CMB angular power spectra for a cosmological model with baryons and an 1 1 eV sterile 
neutrinos (blue line), and for the ACDM model (red line). The points show data from WMAP year 
7 (black), ACT (turquoise) and ACBAR (green). Reproduced from (5]. 



On the contrary, in MOND this massive sterile neutrino would be adequate to 
solve the mass discrepancy on the cluster scale. In fact, the study of the hydro- 
static equilibrium configuration of 30 groups and clusters of galaxies, including the 
two clumps forming the Bullet cluster system, by the analysis of the profiles of the 
galaxy velocity dispersion and hot gas temperature, shows that the mass density pro- 
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file of DM made of 1 1 e V neutrinos always reaches the Tremaine-Gunn limit in the 
cluster center JS] . Some of the dynamical mass must be provided, as expected, by the 
central galaxies, but the required amount implies a mass-to-light ratio of 1 .2Mq /Lq 
in the K band, in agreement with stellar population synthesis models. An 1 1 eV 
sterile neutrino is also consistent with the straight arc, originated by a strong lensing 
effect, observed in the cluster A2390 P9l . 

It is remarkable that the 11 eV sterile neutrino required to match the CMB power 
spectrum also solves the completely independent problem of the mass discrepancy 
of MONDian clusters without requiring any additional free parameter or adjustment 
of the model. It is also intriguing that the Tremaine-Gunn limit is always reached in 
the cluster centers, implying that the dynamical properties of clusters are uniquely 
set by the mass of the sterile neutrino. This relation is completely unkown to the 
standard DM paradigm, where the mass of the cold DM particle does not have any 
role in the dynamical properties of clusters. 

It remains to be seen whether gravitational instability in a universe filled of bary- 
onic matter and one species of sterile neutrino with 11 eV mass can form the ob- 
served cosmic structure at the correct pace, although we remind that the ability to 
explain the cluster mass discrepancy does not directly imply that MOND combined 
with 1 1 eV sterile neutrinos can form clusters in a cosmological context. 

The investigation of structure formation in this model requires an efficient A^- 
body code and proper cosmological initial conditions. Previous A^-body simulations 
of MOND do not consider the presence of any non-baryonic DM and the adopted 
initial conditions might not be consistent for a universe filled with baryons alone. 
These previous A^-body simulations show an overproduction of cosmic structure 
in a baryon-only MONDian universe, in obvious disagreement with observations 

A first attempt of A^-body simulation of structure formation in a MONDian uni- 
verse with baryons and an 1 1 e V sterile neutrino has been performed with a particle- 
mesh cosmological A^-body code that solves the modified Poisson equation of the 
quasi-linear formulation of MOND (QUMOND) and with initial conditions ap- 
propriate to DM made of 11 eV neutrinos 0. The simulation evolved a box of 
5\2h^^ Mpc on a side with 256^ particles from redshift z = 250 to the present time. 

This MONDian hot DM model does indeed produce galaxy clusters with the 
correct order of magnitude of the abundance of observed X-ray clusters, unlike hot 
DM models with standard gravity that have structure formation suppressed on small 
scales ||29l . However the model overproduces the X-ray luminous clusters by a fac- 
tor of three and underproduces the low luminousity clusters with T < 4.5 keV by 
at least a factor of ten. Nevertheless the density profiles of the simulated clusters 
are compatible with the observed profiles of MONDian clusters. In addition, the 
frequency of relative velocities larger than 3000 km s^' of cluster pairs is large 
enough to make likely, unlike the standard model lISOl . the occurrence of systems 
like the Bullet cluster 

Overall the results of this simulation are somewhat unsatisfactory. However, the 
simulation has the very poor mass and length resolutions ~ lO'^M© and ^ 2 Mpc, 
respectively; therefore clustering on small scales is artificially suppressed and this 
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can partly explain the severe underestimate of the abundance of low-luminous X- 
ray clusters. The underproduction of low massive clusters could also be quelled 
by swapping the 11 eV sterile neutrino for a more massive, up to 1 keV, sterile 
neutrinos, because the free-streaming scale decreases with increasing mass of the 
DM particle. As we mentioned above, however, these higher mass sterile neutrinos 
would not be fully thermalised prior to decoupling. The overproduction of massive 
clusters and supercluster size objects actually is a problem that MOND has with any 
mass of sterile neutrino or any other DM particle. 

In addition, the expansion of the Universe is assumed to coincide with the stan- 
dard Friedmann-Robertson-Walker model with a cosmological constant. It is doubt- 
ful that a self-consistent relativistic version of MOND will be close enough to the 
standard model to be consistent with current observational limits on the expansion 
history of the Universe, but discrepant enough to yield the correct structure forma- 
tion rate. Structure formation is affected by the acceleration scale oq which sets the 
enhancement of the intensity of the gravitational field compared to standard gravity. 
Currently, we do not have any observational constraints on the redshift dependence 
of flo, but if ao is lower at higher redshift, the overproduction of cosmic structure 
could be suppressed and the abundance of massive clusters could become consis- 
tent with observations. However, this feature would suppress galaxy formation and 
is therefore not ideal. 

Clearly, the phenomenological formulation of MOND is not conceived to yield a 
self-consistent cosmological model, and a covariant model including MOND needs 
to be implemented to test its predictions of structure formation robustly. For exam- 
ple, the linear perturbation theory of the density field in TeVeS has already been 
outlined 11321 . but the non-linear evolution of cosmic structure, with or without a 
non-baryonic DM component, still needs to be investigated. 



6 Conclusion 

MOND is based on the ansatz that Newonian dynamics is modified when the grav- 
itational field determined by the distribution of baryonic matter drops below the 
acceleration oq ~ 10^'" m s^^ 1941 . The observational evidence for the existence 
of this acceleration scale has been rapidly accumulating over the last decade. The 
success of MOND at describing the dynamics of self-gravitating systems up to the 
scale of galaxies with fewer parameters than the standard DM paradigm is undeni- 
able: in the standard model, the additional gravitational force that is required by the 
observed kinematics is supplied by the presence of proper amounts of DM, but this 
solution requires, unlike MOND, a number of fine tunings and coincidences to ex- 
plain the existence of the acceleration scale that emerges at different length scales. 
MOND provides a far simpler explanation of these observations. 

In addition, on the scale of galaxies, MOND has an impressive predictive power 
that is alien to the DM paradigm. Numerous observations predicted by MOND were 
confirmed years later One of the most striking was that, if interpreted in Newtonian 
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dynamics, the kinematics of LSB galaxies imply that DM dominates the dynamics 
of these systems more than in any other system in the Universe |[95l . 

However, on the scale of galaxy clusters, MOND still requires some form of DM, 
although not as much as in the DM paradigm |7i|. Reproducing the power spectrum 
of the CMB beyond the second peak also requires some non-baryonic DM 1881 . 

In principle, neutrinos might represent a possible candidate for this non-baryonic 
DM. However, current upper limits on the ordinary neutrino mass show that they 
would not be massive enough to condense in the core of clusters and preserve the 
height of the third peak of the CMB power spectrum, as required by observations. 

A viable alternative candidate is a sterile neutrino. This particle seems to be re- 
quired to interpret a number of anomalous results from neutrino experiments ll j. A 
sterile neutrino with mass in the range 11 eV — 1 keV would ensure DM hot enough 
to stream out of galaxies but cold enough to cluster on the scale of massive galax- 
ies and beyond. It is thus in principle capable of explaining the phenomenology of 
astrophysical systems preserving the success of MOND on small scales and, in prin- 
ciple, the success of the standard model on large scales [3|. Such a model, where 
we have a modification of gravity and an exotic DM particle, might be more attrac- 
tive than the current standard ACDM model, because of its elegance on the scale 
of galaxies, the physical motivation of the existence of its DM particle from exper- 
iments on Earth and, eventually, the fact that it requires fewer free parameters than 
ACDM. However, it remains to be investigated whether this model can reproduce 
the full phenomenology of the large-scale structure formation and evolution. 

Moreover, this model rests on MOND, that is still a classical theory, not a covari- 
ant theory. The parent covariant theory that gives MOND in the proper limit is still 
unknown. Many possibilities have been proposed, but none of them appears yet to 
be fully convincing. 

The exciting conclusion is that, whatever model is going to be correct, the exis- 
tence of the acceleration scale, that the data on galaxy scales indicates with astonish- 
ing regularity, has to be explained naturally by the correct model, unless a number 
of surprising coincidences happen to fool us. 
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